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ABSTRACT 

Superbursts are very energetic Type I X-ray bursts discovered in recent years by long term monitoring of X- 
ray bursters, believed to be due to unstable ignition of carbon in the deep ocean of the neutron star. A number of 
intermediate duration bursts have also been observed, probably associated with ignition of a thick helium layer. 
We investigate the sensitivity of these long X-ray bursts to the thermal profile of the neutron star crust and core. 
We first compare cooling models of superburst lightcurves with observations, and derive constraints on the 
ignition mass and energy release. Despite the large uncertainties associated with the distance to each source, 
these parameters are quite well constrained in our fits. For the observed superbursts, we find ignition column 
depths in the range 0.5-3 x 10 12 g cm" 2 , and energy release «2x 10 17 erg g" 1 . This energy release implies 
carbon fractions of X c > 10%, constraining models of rp-process hydrogen burning. We then calculate ignition 
models for superbursts and pure helium bursts, and compare to observations. We show that achieving unstable 
ignition of carbon at accretion rates less than 0.3 of the Eddington rate requires Xc > 0.2, consistent with our 
lightcurve fits. Most importantly, we find that when Cooper pairing neutrino emission in the crust is included, 
the crust temperature is too low to support unstable carbon ignition at column depths of ~ 10 12 g cm -2 . Some 
additional heating mechanism is required in the accumulating fuel layer to explain the observed properties of 
superbursts. If Cooper pair emission is less efficient than currently thought, the observed ignition depths for 
superbursts imply that the crust is a poor conductor, and the core neutrino emission is not more efficient than 
modified URC A. The observed properties of helium bursts support these conclusions, requiring inefficient crust 
conductivity and core neutrino emission. 

Subject headings: accretion, accretion disks-X-rays:bursts-stars:neutron 



1. INTRODUCTION 

In the past few years a new regime of nuclear burning on 
the surfaces of accreting neutron stars has been revealed by 
the discovery of X-ray superbursts (Cornelisse et al. 2000; 
Strohmayer & Brown 2002; Kuulkers 2003). These are rare 
(recurrence times « 1 year), extremely energetic (energies 
« 10 42 ergs) and long duration (4-14 hours) Type I X-ray 
bursts that have been discovered with long term monitor- 
ing campaigns by BeppoSAX and RXTE. Whereas normal 
Type I X-ray bursts involve unstable thermonuclear ignition 
of hydrogen and helium (see Lewin, van Paradijs, & Taam 
1993, 1995; Strohmayer & Bildsten 2003 for reviews), su- 
perbursts are thought to involve ignition of carbon in a much 
thicker layer (Cumming & Bildsten 2001, hereafter CB01; 
Strohmayer & Brown 2002). 

Theoretical studies of superbursts initially focused on their 
potential as probes of nuclear physics. The fuel for su- 
perbursts is thought to be produced by hydrogen and he- 
lium burning at lower densities by the rp-process (Wallace & 
Woosley 1981), a series of proton captures and beta-decays 
on heavy nuclei close to the proton drip line. This pro- 
cess naturally explains the ks 100 s extended tails observed 
in some X-ray bursts (e.g. from the regular burster GS 1826- 
24; Galloway et al. 2004). The amount of carbon remaining 
after H/He burning depends on the details of the rp-process 
(Schatz et al. 2003b), which involves unstable heavy nuclei 
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whose properties are not well known experimentally (Schatz 
et al. 1998). 

CB01 argued that the heavy elements made by the rp- 
process will make the accumulating layer less conductive to 
heat, increasing the temperature gradient within it, and lead- 
ing to earlier ignition than pure carbon models (as had been 
considered earlier by Woosley & Taam 1976, Taam & Pick- 
lum 1978, Lamb & Lamb 1978, and Brown & Bildsten 1998), 
in better agreement with observed superburst energies. In 
fact, Brown (2004) and Cooper & Narayan (2005) showed 
that the ignition conditions are much more sensitive to the 
thermal properties of the neutron star interior, specifically the 
neutrino emissivity of the neutron star core and composition 
of the crust. This is exciting because it offers a new way to 
probe the neutron star interior, complementary to observations 
of transiently-accreting neutron stars in quiescence (Brown et 
al. 1998, Colpi et al. 2001, Rutledge et al. 2002; Wijnands et 
al. 2002; Yakovlev et al. 2004), or cooling isolated neutron 
stars (see Yakovlev & Pethick 2004 and Page et al. 2005 for 
recent reviews). 

Several other long duration X-ray bursts have been ob- 
served that are intermediate in duration and energy between 
normal Type I X-ray bursts and superbursts (e.g. Figure 1 
of Kuulkers 2003). These intermediate bursts have durations 
of w 30 mins and energies w 10 41 ergs, and sources include 
SLX 1737-282 (in 't Zand et al. 2002), 1RXS J171824.2- 
402934 (Kaptein et al. 2000), and 2S 0918-549 (in 't Zand 
et al. 2005). Long duration bursts are expected from accre- 
tion of hydrogen and helium at low rates w 0.01 Mndd (Fu- 
jimoto, Hanawa, & Miyaji 1981; Bildsten 1998; Narayan & 
Heyl 2003; Cumming 2003b). For accretion of solar com- 
position material at these accretion rates, a massive layer of 
pure helium accumulates and ignites beneath a steady hydro- 
gen burning shell. Whereas the hydrogen shell is heated by 



2 



hot CNO hydrogen burning, the helium shell is heated mainly 
by the heat flux emerging from the crust. Therefore, just like 
superbursts, these bursts are potentially sensitive to the crust 
composition and core temperature. The case of pure helium 
accretion is particularly interesting because heating by hydro- 
gen burning then plays no role, making the ignition conditions 
directly sensitive to interior physics. 

It is therefore important to constrain the ignition depth, re- 
currence times, and energy released during superbursts and 
other long X-ray bursts. Our knowledge of superburst recur- 
rence times is limited because they are rare events. Three 
superbursts were seen from 4U 1636-54 separated by in- 
tervals of 2.9 and 1.8 years (Wijnands 2001; Kuulkers et 
al. 2004). Dividing the total duration of observations of X- 
ray bursters with the BeppoSAX/WFC by the number of su- 
perbursts observed gives a recurrence time estimate of 0.4-2 
years (in 't Zand et al. 2003). Brown (2004) and Cooper & 
Narayan (2005) emphasized that to achieve ignition of carbon 
on w 1 year timescales at accretion rates M w 0.1 M^dd re- 
quires the accumulating layer to be sufficiently hot. An en- 
hanced core neutrino emissivity (as would be produced if, 
for example, the direct URCA process operated in the core, 
e.g. Yakovlev & Pethick 2004, Page et al. 2005) together with 
a large crust conductivity gives very long (;§> 10 yr) superburst 
recurrence times, inconsistent with observations. Recently, 
in 't Zand et al. (2005) showed that the intermediate dura- 
tion X-ray burst from 2S 0918-549 is well explained by ac- 
cretion of pure helium at the observed rate of M w 0.01 Mem, 
assuming that most of the heat released in the crust by pyc- 
nonuclear reactions and electron captures (Haensel & Zdunik 
1990, 2003) flows outwards and heats the accumulating he- 
lium layer. However, they did not explore the sensitivity of 
this assumption to the interior physics. 

In this paper, we investigate the constraints on interior 
physics coming from superbursts and pure helium bursts. We 
first derive independent constraints on the ignition depth and 
energetics of superbursts by fitting the observed lightcurves to 
theoretical cooling models as calculated by Cumming & Mac- 
beth (2004, hereafter CM04). We then calculate ignition con- 
ditions for both superbursts and pure helium bursts, and com- 
pare with observed properties. We start in ^2]by summarizing 
the properties of our cooling models for superbursts, present 
the fits to the observed lightcurves, and discuss the constraints 
on the ignition depth and energy release. In §0 we calculate 
ignition conditions for superbursts and discuss the implica- 
tions for the thermal structure of the interior. We show that the 
best fit is obtained for inefficient neutrino emission in the neu- 
tron star crust and core. Next, in |4] we apply these ignition 
models to pure helium bursts, and show that their properties 
imply the same conclusion: inefficient neutrino emission. We 
conclude in |5] In Appendix A, we discuss a simple model of 
the early phase of the superburst lightcurve which reproduces 
the time-dependent results, and in Appendix B give approxi- 
mate analytic solutions for the crust temperature profile. 

2. COOLING MODELS FOR SUPERBURSTS AND 
COMPARISON TO OBSERVATIONS 

2.1. Cooling models 

CM04 computed lightcurves for superbursts by assuming 
that the fuel is burned locally and instantly at each depth, and 
then following the thermal evolution and surface luminosity 
as the burning layers cool. They showed that the lightcurve 
of the cooling tail of the superburst is a broken power law, 
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Column depth (g cm -2 ) 

FIG. 1 . — Energy radiated from the surface after 3 hours (dotted curves), 6 
hours (solid curves), and 12 hours (dashed curves) for £"17 = 1, 1.5,2,2.5 and 
3, as a function of the column depth. We assume a neutron star radius R = 10 
km(£rad.oo ocR 2 ). 



with time of the break proportional to the cooling time of the 
entire layer. The early phase of cooling depends mostly on 
the energy released in the flash; the late phase of cooling de- 
pends mostly on the thickness of the layer. We now apply 
these models to the observations, and discuss the constraints 
on superburst ignition conditions. 

We refer the reader to CM04 for full details of the cooling 
models, including approximate analytic expressions for the 
flux as a function of time. The parameters of the models are 
the ignition column depth (y^ 10 12 g cm" 2 ), and the energy 
release per gram (En 10 17 ergg" 1 ) which is assumed to be 
independent of depth. The thermal evolution is followed nu- 
merically by the method of lines, which involves finite differ- 
encing on a spatial grid, and then integrating the resulting set 
of coupled ordinary differential equations for the temperature 
at each grid point forward in time. It is important to note that 
the rise of the superburst is not resolved, since the entire fuel 
layer is assumed to burn instantly. The models assume a neu- 
tron star mass and radius of M = 1 .4 M Q and R = 1 km, giving 
the surface gravity g = (GM/R 2 )(\+z) = 2.45 x 10 14 cm s" 2 
and redshift 1 + z = 1.31. We include the equation of state, ra- 
diative and conductive opacities, heat capacity, and neutrino 
emissivities as described by Schatz et al. (2003a). At the tem- 
peratures and densities appropriate for a superburst, the neu- 
trino emission is mostly due to pair annhilation (CM04). 

We first summarize some of the properties of the models. 
CM04 already noted that the power law decay gives a long 
tail to the superburst lightcurves, which is similar to the long 
tails observed in superbursts (Kuulkers et al. 2002; Cornelisse 
et al. 2002). Figure Q shows the amount of energy radiated 
from the surface in the first 3, 6, and 12 hours as a function of 
column depth for different choices of £17. The insensitivity of 
radiated energy to column depth fory > 10 12 g cm" 2 in Figure 
^ shows that the total emitted energy is not a good indicator 
of the ignition column depth. There are two reasons for the 
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FIG. 2. — The "thermostats" of neutrino emission and inwards conduction 
of heat. Upper panel: ratio of energy released as neutrinos to energy radiated 
from the surface, both in the first 24 hours. Neutrinos dominate the energy 
release for large y and £17. Lower panel: Fraction of the total nuclear energy 
released that escapes in the first 24 hours, either as neutrinos or from the 
surface. For large y, a significant fraction of the energy released is conducted 
inwards and released on a longer timescale. 



characteristic radiated energy of w 10 42 ergs (Strohmayer & 
Brown 2002). First, neutrino emission takes away most of the 
energy for large columns, and secondly heat flows inwards to 
be released on longer timescales. These effects are quantified 
as a function of y\2 and En in Figure|3] The first panel shows 
the ratio of energy lost as neutrinos to the energy lost through 
the surface. For example, for y k, 10 13 g cm" 2 and En k, 3, 
neutrinos take away an order of magnitude more energy than 
is lost from the surface. This is in rough agreement with the 
one-zone model of Strohmayer & Brown (2002). The second 
panel shows the fraction of energy that is lost in the first 24 
hours, either as neutrinos, or from the surface. The remaining 
energy, which is released on longer timescales, can be a sig- 
nificant fraction of the total for column depths w 10 13 g cm" 2 



and £17 w 1. 

Even without detailed fits to observed lightcurves, these 
results give some indication of the values of £17 needed 
to match the observed properties of superbursts. Figure ^ 
shows that an energy release £17 > 1 is required for the ob- 
served burst energy to reach > 3 x 10 41 ergs during the first 
few hours. On the other hand, for large values of £17, the 
initial flux exceeds the Eddington flux F-Edd = cg/n = 2.2 x 
10 25 erg cm" 2 s" 1 (gi 4 /2.45)(1.7 /(1 + X)), in which case the 
superburst would be expected to show photospheric radius ex- 
pansion. In Figure we plot the time for which the flux ex- 
ceeds ^Edd for different £17 values. This time is not very sen- 
sitive to the ignition column, since the early evolution of the 
burst is independent of the layer thickness. For £17 > 2, the 
flux is super-Eddington for timescales of minutes or longer. 
The superburst from 4U 1820-30 showed an extended period 
of photospheric radius expansion lasting for ~ 1000 s (SB02). 
This is in good agreement with the expectation that this source 
had a significant energy release due to large amounts of car- 
bon produced by stable burning of pure helium (SB02; Cum- 
ming 2003a). Figure [3] implies that £17 > 5 is required to 
get such a long period of super-Eddington luminosity with 
pure helium. However, there is no strong evidence for pho- 
tospheric radius expansion in any other superburst 5 . Taken 
together, these two constraints imply that £17 « 2 for most 
superbursts. 



2.2. Fits to superburst lightcurves 

We have fitted the superburst lightcurves to the cooling 
models. The parameters of the models are £17 and yn- How- 
ever, there are two additional parameters in our fits. The rise 
of most superbursts is not observed because of data gaps, 
making the start time of the burst uncertain, and so we in- 
clude the start time as an extra parameter. Most importantly, 
the distance to the source is not well constrained in most 
cases. This dominates the uncertainty in the fitted parame- 
ters, and so we have fitted the models by holding distance 
fixed at different trial values, and searching over the remain- 
ing parameters to find the best fitting model at each distance. 
We include BeppoSAX/WFC data for the superbursts from 
4U 1254-690 (in 't Zand et al. 2003), KS 1731-260 (Kuulk- 
ers et al. 2002), 4U 1735-444 (Cornelisse et al. 2000), Ser 
X-l (Cornelisse et al. 2002), GX 17+2 (in't Zand et al. 2004), 
and the RXTE/PCA lightcurve of 4U 1636-54 (Strohmayer & 
Markwardt 2002; Kuulkers et al. 2004). For GX 17+2, we 
use burst 2 from Figure 7 of in 't Zand et al. (2004). This is 
one of the best candidates for a superburst, and has the most 
complete lightcurve. 

We have extended the CM04 models to a large grid in 
£17 and yi2 for comparison to the observations. For a given 
source distance, we calculate the flux at the surface of the 
star F+ which corresponds to the observed peak flux / pea k, 
i.e. 4irR 2 F+ = 47raf 2 / pea k. We will refer to £* in units of 
10 24 erg cm" 2 s" 1 as £24 . This quantity sets the normaliza- 
tion scale for comparison with the theoretical models, and is 



5 Precursors were seen with BeppoSAX/WFC from KS 1731-260, 
4U 1254-69, and GX 17+2. In GX 17+2, the spectral data are of insufficient 
quality to see radius expansion during the precursor or the minutes thereafter 
(because of the high persistent emission); in KS 1731-260, no radius expan- 
sion was seen (Kuulkers et al. 2002); in 4U 1254-69 there are no indications. 
In all cases the peak flux of the precursor is smaller by a factor of 1 .5-2 than 
the brightest of the ordinary bursts. 
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FIG. 3. — Time for which the flux exceeds the Eddington flux as a function 
of energy release £17. The curves are for y = 10 11 (long-dashed), 3 X 10 11 
(short-dashed), 10 12 (dotted), and 10 13 g cm -2 (solid). We show two sets 



of curves for solar composition (/"Edd = 2.2 X 
helium (£ Edd = 3.7 X 10 25 erg cm" 2 s" 1 )- 
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TABLE 1 
Fits to superburstlightcurves 



Source 


/peak 


d/R b 


£n c 


yn c 


4U 1254-690 


0.22 


13 


1.5 


2.7 


4U 1735-444 


1.5 


8 


2.6 


1.3 


KS 1731-260 


2.4 


4.5 


1.9 


1.0 


GX 17+2 burst 2 


0.8 


8 


1.8 


0.64 


SerX-1 


1.9 


6 


2.3 


0.55 


4U 1636-54 


2.4 


5.9 


2.6 


0.48 



"Observed peak flux in units of 10 8 erg cm 2 s '. 
b Adopted distance in units of kpc/10 km. 

c The fitted parameters scale roughly as £17 oc (rf/R) 8 / 7 andyi2 <x (d/R) 10 / 1 
(see text). For a 50% distance uncertainty, the uncertainties in £17 and V12 
are 60% and 70% respectively (see also Fig. [4}. 



given by 



F 24 = 9.5 



,/peak 



10 



" 8 erg cm 2 s 



d/R 



10 kpc/10 km 



(1) 



We then search for the minimum value of \ 2 over the grid of 
theoretical models, with En ranging from 0.5 to 3 in steps of 
0.125, and V12 ranging from 10 11 to 3 x 10 13 g cm" 2 in steps 
of 1/16 in log 10 y (i.e. factors of 15% between successive y 
values). For each model, we vary the start time of the super- 
burst between the limits allowed by the observations to find 
the best fit. Before comparing the model to the data, we red- 
shift the time and flux assuming a gravitational redshift factor 
of 1 +z= 1.31, appropriate for a 10 km, 1.4 M Q neutron star. 
Our results are not very sensitive to variations in the redshift 
factor within the expected range (roughly 1 .2-1 .5). 

Figure 0] shows the best fitting En, y, and the reduced x 2 
of the fit for each source, as a function of both F24 and the 
distance to radius ratio d/R. A larger flux normalization for 
the observed lightcurve results in larger values of En and 
yi2, which increase in such a way as to maintain the over- 



all shape of the cooling curve. The scalings are straightfor- 
ward to understand from the analytic expressions for the flux 
given by CM04 (see eq. [4] of that paper). At early times, 

the flux is F oc t~ - 2 E 7 ^ 4 (independent of column depth). In 
Appendix A, we discuss the physics underlying these scal- 
ings. Comparing with the fitted values to set the constant, we 

4/7 

find £17 » 0. 8F 2 4 > which is in good agreement with the ob- 
served relation between the fitted £17 value and F24. There 
is a similar scaling for the best fit column depth, which can 
also be understood from the analytic fit, but now at late times, 

where F oc yE 1 ^ 2 , giving y oc F 5/ " 7 . For a given fit, the frac- 
tional uncertainties in En and y can therefore be estimated as 
w (A/1)(8F/F) and (5/7)(5F/F) assuming that the distance 
uncertainty dominates. 

For specific choices of distance to each source, we show the 
best fitting models in Figures 5 to 10, and list the parameters 
in Tabled The fitted values can be rescaled to a different dis- 
tance using the analytic scalings, or by referring to Figure 4. 
We adopt distance estimates from the literature for 4U 1254- 
690 (in 't Zand et al. 2003) and 4U 1636-54 (Augustein et 
al. 1998). For GX 17+2 and 4U 1735-444 we adopt a fiducial 
value of8kpc. ForSerX-1 andKS 1731-260, we take a lower 
distance then the upper limits or estimates in the literature, be- 
cause this significantly improves the fit of our models. For ex- 
ample, Muno et al. (2000) place a distance limit of d < 7 kpc 
for KS 1731 -260 using radius expansion X-ray bursts, assum- 
ing that the peak luminosity is the Eddington luminosity for 
pure helium. We find that for d/R>5 kpc/10 km the super- 
burst lightcurve is not well fit by our models, with x 2 rapidly 
increasing for larger d/R. If the distance is 7 kpc, the required 
neutron star radius is > 13 km. Alternatively, the source could 
be closer. For example, using the Eddington luminosity for a 
solar composition rather than pure helium gives a closer dis- 
tance by a factor of (1.7) 1 / 2 or 1.3. We choose a distance 
d/R = 4.5 kpc/10 km for the fit shown in Figure|5] 

The most detailed lightcurve is for4U 1636-54. This source 
has shown three superbursts (Wijnands 2001; Strohmayer & 
Markwardt 2002; Kuulkers et al. 2004), but we show here 
the superburst observed by RXTE/PCA (Strohmayer & Mark- 
wardt 2002). Figure[6]shows the Standard 1 mode lightcurve, 
which has a time resolution of 1 second, but has been binned 
to 10 second resolution for clarity. The best fit model agrees 
well with the observed decay. However, there are differences 
at the w 10% level between the model and the shape of the ob- 
served lightcurve. The BeppoSAX data for the other sources 
have a much lower time resolution, but still allow a good 
constraint on the ignition depth. The fitted column depth 
is not very sensitive to the assumed start time of the super- 
burst. Most important is how quickly the luminosity decays 
away from the peak value. For example, the count rate for 
4U 1254-690 takes several hours to fall to 30% of the peak 
value, whereas for 4U 1636-54 the count rate reaches 30% of 
the peak after less than 2 hours. We fit only to data points that 
have count rates more than 10% of the peak value. At low 
luminosities, the lightcurve is sensitive to how well the accre- 
tion luminosity has been subtracted, however, this uncertainty 
has only a small effect on the fitted column depth. 

The best-fitting models have En in the range 1.5 to 2.6. 
As we argued in §2.1, lower values of En give a luminos- 
ity at early times that is smaller than observed. At the upper 
end of this range, Figure [3] shows that the flux should exceed 
the Eddington flux for several minutes, inconsistent with the 
lack of observed photospheric radius expansion. This may 
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FIG. 4. — Left panel: best fitting En and y, and the associated reduced \ > as a function of assumed peak flux F24- The fitted values approximately follow the 

scalings E„ ss 0.8F,, 4 / 7 and y oc F 2 5 / 7 . We show results for 4U 1254-690 (short-dashed), KS 1731-260 (long dashed-short dashed), 4U 1735-444 (solid), Ser X-l 
(long-dashed), GX 17+2 (burst 2 dot-dashed, burst 3 long-dot-dashed), and 4U 1636-54 (dotted). Right panel: same as left panel, but now using the observed 
peak flux to plot everything in terms of the distance to the source. The \ 2 f° r 4U 1636-54 (dotted curves) is off scale in the lower panel. 
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FIG. 5. — Fitted lightcurve for KS 1731-260, assuming the distance given 
in Table 1. Solid data points are included in the fit, open data points (with 
fluxes less than 0. 1 of the peak flux) are not included. 
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FIG. 6.— Fitted lightcurve for4U 1636-54. 
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indicate that the burning does not extend all the way to the 
surface, which our models assume, but instead stalls at a loca- 
tion where the thermal time to the surface is of order minutes. 
More generally, our models are not valid for times less than 
the superburst rise time. Also, we have not fitted our mod- 
els to the superburst from 4U 1820-30, which was observed 



by RXTE/PCA (Strohmayer & Brown 2002). This superburst 
had a complex lightcurve, with an extended period of photo- 
spheric radius expansion, lasting about 1000 seconds, indicat- 
ing a large energy release. More detailed ID models which 
can follow the superburst rise are needed to address both of 
these issues. 

The best-fitting column depths are in the range 0.5-3 x 
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FIG. 7.— Fitted lightcurve for4U 1254-690. 
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FIG. 8.— Fitted lightcurve for4U 1735-444. 



10 12 g cm" 2 . Larger column depths closer to 10 13 g cm" 2 are 
not consistent with the observed lightcurves. We have ranked 
the sources in Tabled in order of decreasing column depths. 
This ordering approximately reproduces the ordering of su- 
perbursts by their observed durations in Figure 7 of in 't Zand 
et al. (2004). As pointed out in that paper, the superbursts 
from the rapidly accreting source GX 17+2 have low column 
depths, but not significantly lower than other superbursts. We 
find that the GX 17+2 burst has a similar ignition depth to the 
superbursts seen from 4U 1636-54 and Ser X-l. 

CM04 derived constraints on the ignition column depth 
from the observed "quenching" of normal Type I bursting 
behavior for weeks following a superburst (e.g. Kuulkers et 
al. 2002). The layer continues to cool well after the superburst 
luminosity falls below the accretion luminosity. This residual 
heat flux quenches the instability of H/He burning (CB01). 




Time (hours) 



FIG . 9 . — Fitted lightcurve for Ser X- 1 . 




Time (hours) 



FIG. 10.— Fitted lightcurve for GX 17+2 (burst 2 from in 't Zand et 
al. 2004). 



CM04 showed that the observed limits on the quenching 
timescale, although not very constraining, were at least con- 
sistent with ignition at column depths of « 10 12 — 10 13 g cm" 2 . 
Our fitted values of yn are consistent with Figure 4 of CM04, 
except for KS 1731-260. The quenching timescale implies 
yn 3 for this source, whereas our fit gives yi2 ~ 1. Apart 
from the uncertainties associated with distance, one possibil- 
ity is that the flux required to stabilize H/He burning is a factor 
of 3 lower than the crude estimate of CM04 (their eq. [5]). 

3. IGNITION MODELS FOR SUPERBURSTS 

The fits to the superburst lightcurves imply ignition column 
depths w (0.5-3) x 10 12 g cm" 2 , which is accumulated in 2-10 
years at 0.1 rhm&, or 0.6-3 years at 0.3 WEdd, roughly consis- 
tent with the observational constraints on recurrence times. In 
this section, we compare our fits to ignition models for su- 
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perbursts, and use the cooling models to predict superburst 
properties as a function of accretion rate. In particular, we 
emphasize the constraints on the thermal profile of the crust 
and temperature of the core. 

3.1. Details of the ignition calculations and physics input 

Our ignition models follow those of Brown (2004), and are 
extensions of the CB01 carbon ignition models. However, 
we now integrate down to the crust/core interface, solving the 
thermal structure of the crust directly, rather than taking the 
outwards flux from the crust as a free parameter. Following 
the simplifications of Yakovlev & Haensel (2003) and Brown 
(2004), we do not integrate the full structure of the star, but 
adopt a plane-parallel approximation, and take the gravity g 
and gravitational redshift factor 1 +z to be constant across the 
crust. Our independent variable is then the column depth y, 
where hydrostatic balance gives the relation y = P/g (units 
of mass per unit area). We integrate the heat equation and 
entropy equation 

d 2L— 4. n\ 

, — ^nuc ' £v \£) 

dy 



dT _ F 
~dy~~pK 



(3) 



where K is the thermal conductivity, F the heat flux, T the 
temperature, e v the neutrino emissivity, and e nuc the rate of 
heating from nuclear reactions in the crust. All quantities in 
equations and (0 are proper (local) quantities if we inter- 
pret y as the rest mass column depth (e.g. see §4.1 of Cum- 
ming et al. 2002). 

The flux flowing outwards from the crust F out is a crucial pa- 
rameter, since this flux heats the accumulating fuel layer. We 
write it in terms of the parameter Qb as F out = rnQb- Follow- 
ing Brown (2000), we model the heating from pycnonuclear 
reactions and electron captures in the crust by setting e nuc = 
Qnu C rn/Ay to be constant over the range Ay from y = 5.7 x 
10 15 to 2.2 x 10 17 g cm" 2 , where Q mc = 1 .4 MeV per nucleon 
is the total energy release in the crust. 

In the crust, the equation of state is determined by degen- 
erate electrons and neutrons. We use the results of Mackie & 
Baym (1977) to correct the neutron Fermi energy for inter- 
actions. The top of the crust is set by T = Z 2 e 2 jak^T = 175 
(Potekhin & Chabrier 2000), where a = (3/47m,-) 1/3 is the in- 
terion spacing. Since the pressure is dominated by relativistic, 
degenerate electrons, this translates to a column depth 



ymeit = 4.5 x 10 13 gcm 



aV /3 

56 ) 



(4) 



where 7g = T/10 8 K. For the core, we adopt the sim- 
plified equation of state p{r) = p c (l - (r/R) 2 ), where M = 
(87r/15)(p e ^ 3 ) (Yakovlev & Haensel 2003). This gives p c = 
1.7 x 10 15 g cm" 3 (M/1.4 M )(/?/lO km)" 3 . For a density 
at the crust/core interface of 1.6 x 10 14 gem" 3 (Lorentz et 
al. 1993), the core radius is w 0.95R. In fact, in our mod- 
els, the depth of the core/crust interface from the surface is 
w 1 km w 0.1/?. We expect this difference to have only a 
small effect on our results. 

We integrate the temperature profile inwards, changing the 
composition from "fuel" to "ash" at a depth y lga . We set the 
outer boundary at y = 10 8 g cm" 2 , and take the temperature 
there to be 2 x 10 8 K. We have checked that the ignition depth 



is not very sensitive to this choice of outer temperature 6 . We 
iterate to find the choice of flux at the surface that results in 
the flux at the base matching the core neutrino luminosity, 
F c +L u (T c )/AitR 2 = 0. We write the core luminosity as roughly 
L v » (A-k /J)R i Q u , where Q v is the emissivity per unit vol- 
ume, giving the inner boundary condition F c = —RQ„(T c )/3. 

We calculate the ignition criterion for carbon according to a 
one-zone approximation (Fujimoto, Hanawa, & Miyaji 1981; 
Fushiki & Lamb 1987b; Cumming & Bildsten 2000), compar- 
ing the temperature sensitivity of the heating rate at the base 
of the layer to the that of a local approximation to the cool- 
ing rate. Note that we calculate the temperature sensitivity of 
the heating rate numerically rather than assume a particular 
value (Brown 2004 assumed that d In ecc/d\n T = 26). Al- 
though approximate, this ignition criterion agrees well with 
more detailed linear stability (Narayan & Heyl 2003) and 
time-dependent calculations (Woosley et al. 2004) for H/He 
burning, and we expect it to be accurate here also (Cooper & 
Narayan 2005). The carbon burning rate is given by Caughlan 
& Fowler (1988) with screening from Ogata et al. (1993) 7 . In 
addition, at the ignition point we check whether the timescale 
for carbon depletion is longer than the accumulation time. As 
shown by CB01, carbon burns stably during accumulation for 
low accretion rates, so that the thermal instability is avoided. 
We show only results for which the depletion time is longer 
than the recurrence time. Our results compare well with those 
of Brown (2004). Typically, the recurrence times we find are 
a factor of < 50% larger than those in Brown (2004), after 
correcting those results for gravitational redshifting. 

3.2. Neutrino cooling and crust composition and 
conductivity 

The main parameters in our models are the crust neutrino 
emissivity, the core neutrino emissivity, and the crust compo- 
sition and conductivity. In the crust, we include cooling due 
to neutrino Bremsstrahlung according to Haensel, Kaminker, 
& Yakovlev (1996) in the liquid phase, and Kaminker et 
al. (1999) in the solid phase. The fitting formula given by 
Kaminker et al. (1999) is for an equilibrium crust compo- 
sition. To account for an accreted composition, we multi- 
ply the emissivity by a factor R where log l() R = -0.2 for 
p < 10 11 g cm" 3 , -0.3 for 10 11 g cm" 3 < p < 10 13 g cm" 3 , 
and -0.4 for p > 10° g cm" 3 . This closely reproduces the 
results for accreted matter shown in Figure 7 of Kaminker et 
al. (1999), and agrees to a factor of 3 with the formula of 
Haensel et al. (1996) for densities p > 10 12 g cm" 3 . 

Most importantly, we include the possibility that the neu- 
trons in the crust are superfluid. In this case, there is an addi- 
tional neutrino cooling mechanism involving the continuous 
formation and breaking of Cooper pairs (Flowers, Ruderman, 
& Sutherland 1976; Voskresensky & Senatorov 1987). We 
use the emissivity cal culat ed by Yakovlev, Kaminker, & Lev- 
enfish (1999) (see eq 1C5I ). and we take the neutron l So crit- 
ical temperature T c as a function of density as given by the 
calculation of Schwenk, Friman, & Brown (2003). We have 
also used the results of Ains worth, Wambach, & Pines (1989) 
(case 2 from their Fig. 3), which has a slightly different profile 

6 The most sensitive case is for fast cooling in the core, and Q = 100 in the 
crust, for which changing the outer temperature by a factor of 2 increases the 
outwards flux by a factor of 2, and increases the ignition depth by 5%. 

7 The screening calculations of Ogata et al. (1993) are not appropriate for 
the pycnonuclear regime (e.g. Kitamura 2000; Gasques et al. 2005). How- 
ever, for the temperatures Tg > 2 that we consider in this paper, carbon burn- 
ing is safely thermonuclear. 
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TABLE 2 
Core neutrino emission 



Label 


Type a 


Prefactor b 
(erg citT 3 s" 1 ) 


Comment 


a 


fast 


10 26 


fast cooling 


b 


slow 


3 X 10 21 


enhanced 


c 


slow 


10 20 


mURCA 


d 


slow 


10 19 


nn Bremsstrahlung 


e 


slow 


10 17 


suppressed 



"Fast and slow cooling laws are of the form Q v = Q/(T c /\(f K) 6 and Q v = 
Q s {T c /\(f K) 8 respectively. 
b Either Q s or Qt for slow or fast cooling, respectively. 
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FIG. 1 1 . — The effect of core neutrino emissivity on superburst ignition 
conditions at m = 0.3 tiiEdd- We assume a disordered lattice in the crust, 
and do not include Cooper pairing. The accreted composition is 20% 12 C 
(Xc = 0.2) and 80% 56 Fe by mass. From top to bottom, the temperature 
profiles are for increasing core neutrino emissivity; the letters refer to Table 
|z| The long-dashed line shows the carbon ignition curve for Xc = 0.2, and the 
vertical dotted line indicates a column depth of 10' 2 g cm -2 . 



and a larger maximum temperature, but the results are simi- 
lar and so we do not show them here. Cooper pair emission 
was not considered by Brown (2004) and Cooper & Narayan 
(2005); however we show here that it has a dramatic effect on 
the crust temperature profile. 

For the core neutrino emissivity, we consider the "fast" 
and "slow" cooling laws Q v = Q f (T c /l0 9 K) 6 and Q v = 
Q S (T C /\Q 9 K) 8 (e.g. Yakovlev & Haensel 2003; Yakovlev & 
Pethick 2004, Page et al. 2005). The "standard" slow cool- 
ing by modified URCA processes has Q s ~ 10 20 erg cm" 3 s" 1 . 
However, if either the core protons or neutrons are super- 
fluid, with very high values of T c (3> 10 9 K), then this pro- 
cess is totally suppressed, leading to cooling by nucleon- 
nucleon Bremsstrahlung (involving the non-superfluid com- 
ponent). This process is roughly a factor of ten slower than 
modified URCA, and so we take Q s ~ 10 19 erg cm" 3 s" 1 in 
this case. If both protons and neutrons are strongly super- 
fluid in the core, the neutrino emission will be supressed 
further. To model this case, we assume that the core neu- 
trino emission is suppressed by a further factor of 100, giving 
Q s ~ 10 17 erg cm" 3 s" 1 . However, in the more reasonable case 



\ 


I I I 1 , 


, , , 


i i , i . 


" \ 












^ - 




\ 






_ ~\ d - 


\ 








- \ 








' ^> 


















\ 




\ \ 

\ ^ 




\ 








\ 
























, , i 


, i 


\ \ " 

, , 1 , 



10 12 14 16 18 

log Column Depth (g cm -2 ) 



FIG. 12. — The effect of crust composition and conductivity on superburst 
ignition conditions. Temperature profiles for superburst ignition models at 
m = 0.3 WEdd. We show two cases of core neutrino emissivity: slow cooling 
with Q s = 10" erg cm -3 s~' and fast cooling with Qf = 10 26 erg cm" 3 s" 1 . 
Solid lines are for a composition of 56 Fe and a disordered lattice. Short- 
dashed lines have a heavier composition (A = 106, Z = 46), and dot-dashed 
lines are for a larger thermal conductivity (Q = 100). The long-dashed line 
shows the carbon ignition curve for Xc = 0.2, and the vertical dotted line 
indicates a column depth of 10 12 g cm -2 . 



that the neutron and/or proton T c in the core are of the order 
of 10 9 K there is intense neutrino emission from the Cooper 
pair formation, resulting in an enhanced slow cooling rate 
which we model by considering Q s ~ 3 x 10 21 erg cm" 3 s" 1 
(see, e.g., Figures 20 and 21 in Page et al. 2004). Finally, we 
also consider a fast cooling rate with Qf ~ 10 26 erg cm" 3 s" 1 
corresponding, e.g., to the direct Urea process. These mod- 
els are summarized in Table |2] The core temperature T c 
can be estimated in each case. For slow cooling, we find 

T c « 4.9 x 10 8 K (fl /H /Ql%)(m/m EM y /& and fast cooling 

T c « 5.0 x 10 7 K (fU 6 /Q f %) (m/m E dd) 1/6 where f in is the 
fraction of heat released in the crust that is conducted into the 
core. 

For the composition of the crust, we use the composition 
calculated by either Haensel & Zdunik (1990) or Haensel & 
Zdunik (2003). The difference between these two calcula- 
tions is the nucleus assumed to be present at low densities, ei- 
ther 56 Fe (Haensel & Zdunik 1990), or a heavy nucleus 106 Pd 
(Z = 46) (Haensel & Zdunik 2003), as would be appropriate 
if rp-process hydrogen burning is able to run to its endpoint 
(Schatz et al. 2001). We calculate results for these two cases 
to illustrate the variation expected from changes in composi- 
tion. For the conductivity, we consider two cases. The first 
is a "disordered" crust, for which we take the conductivity 
to be that of a liquid phase, in the second case, we calculate 
the contributions fromphonons (Baiko & Yakovlev 1996) and 
electron-impurity scattering (Itoh & Kohyama 1993), taking 
the impurity parameter Q = 100 (see Itoh & Kohyama 1993 for 
a definition of the impurity parameter, written as ((AZ) 2 ) in 
their notation). Note that a crust with Q = 100 is very impure. 
However, we do not consider smaller values of the impurity 
parameter because as we will show they would not agree with 
observed X-ray burst properties. 
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FIG. 13. — The effect of neutrino cooling by Cooper pairs on superburst 
ignition conditions. For two different core neutrino emissivities, we show 
temperature profiles with (solid) and without (dot-dashed) neutrino cooling 
by Cooper pairs. These models are for m = 0.3 HiEdd an d nave = 0.2. 
The long-dashed line shows the carbon ignition curve for Xq — 0.2, and the 
vertical dotted line indicates a column depth of 1()' 2 g cm -2 . 



3.3. Ignition conditions at a fixed accretion rate 

We first calculate ignition conditions for carbon at m = 
0.3 «iEdd, for different values of neutrino emission and crust 
properties. In order to illustrate the effects of different param- 
eters, we start by assuming that the neutrons in the crust are 
normal (no cooling due to Cooper pair neutrinos). For this 
case, Figure ^2 shows the effect on the temperature profile 
of varying the core neutrino emissivity. Less efficient neu- 
trino emission leads to a higher core temperature and a greater 
fraction of the energy released in the crust is emitted through 
the surface, heating the carbon layer. At this accretion rate, 
standard slow cooling in the core results in recurrence times 
> 3 years, longer than inferred from observations. Some sup- 
pression of the modified URCA rate is necessary to bring the 
predicted and observed recurrence times into agreement. 

Figure^lillustrates the effect of the crust composition and 
conductivity on the temperature profile for fast and slow cool- 
ing in the core. At m = 0.3 rnEdd, the model with Q = 100 and 
fast core neutrino emission gives superburst recurrence times 
that are much longer than observed. In general, the change 
in ignition depth with composition is much smaller than the 
change in ignition depth with other model parameters. This is 
because a heavier composition decreases the thermal conduc- 
tivity, but also results in less outwards flux from the crust, as 
pointed out by Brown (2004) and Cooper & Narayan (2005). 
In Appendix B we give a simple analytic argument to under- 
stand this. 

We now include neutron superfluidity in the crust. Figure 
EJshows the dramatic effect of the extra cooling from Cooper 
pair neutrino emission. We show profiles for either fast cool- 
ing in the core, or highly suppressed core cooling, with and 
without Cooper pair cooling in the crust. When Cooper pair- 
ing of neutrons is included, the temperature is limited by neu- 
trino losses to a value T < 5 X 10 K. This is true even for 
highly suppressed neutrino cooling in the core; in this case, 
most of the energy release in the crust leaves as neutrino emis- 
sion from within the crust itself. In Appendix B, we show 




log Density (g cm 3 ) 

FIG. 14. — The critical temperature for neutron superfluidity in the crust, 
according to Schwenk et al. 2003, and an example of the neutrino emis- 
sivity as a function of depth from Cooper pairing (solid lines) and electron 
brems strah lung (dotted line), for model "e" with Cooper pairing shown in 
Figure [T3l The low density peak in the Cooper pair emissivity corresponds 
to column depths in the range 10 16 -10 17 g cm -2 , the peak at higher densities 
corresponds to column depths 3 X 10 18 g cirT 2 near the base of the crust. 



how to understand this limiting temperature analytically by 
balancing the heating and cooling rates. 

Figure ^] shows that if Cooper pairing is important in 
the crust, the superburst recurrence times should be long 
w 10 years, and insensitive to core neutrino emissivity. In 
Figure ^] we show the critical temperature for the neutron 
superfluid, and the neutrino emissivity from Cooper pairing 
and electron Bremsstrahlung for comparison. The Cooper 
pair emission is concentrated in two regions where T ~ T c , as 
discussed by Yakovlev et al. (1999), and is therefore sensitive 
to the behavior of T c close to the superfluid threshold, which 
is uncertain. We address these uncertainties in Appendix C. 
However, we find that as long as the critical temperature in- 
creases from zero to large values, crossing the crust temper- 
ature, the process is important. For example, we have tried 
modelling the T c profile as log-Gaussian in density, and vary- 
ing the central density and width, but have not been able to 
significantly reduce the Cooper pair neutrino emission. The 
peak in emission at lower densities has the largest effect, since 
this extra cooling occurs at the location of the energy release 
in the crust (close to neutron drip). The peak in emissivity 
near the core boundary has a smaller effect, equivalent to an 
extra core neutrino emission. 

3.4. Variation with accretion rate and comparison to 
observations 

CB01 showed that the ignition conditions are very sensi- 
tive to accretion rate, and so a natural question is how much 
the results of §3.3 depend on the assumed accretion rate. We 
explore this dependence here. The accretion rate in the super- 
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FIG. 15. — Flux from the crust heating the fuel layer and ignition column 
depth as a function of accretion rate. The solid curves show results for a dis- 
ordered crust, and a composition of Xc = 0.2 and 56 Fe, for the four different 
core neutrino emissivities of Figure 1111 More efficient core neutrino emis- 
sion gives a lower flux from the crust Q/,, and a larger ignition column depth. 
Other curves show variations on the "d" model. The short-dashed curve is 
for a heavier composition (A = 106, Z = 46); the dot-dashed curve is a higher 
crust conductivity (Q = 100); the dotted curve includes Cooper pairing in the 
crust (labelled "SF"). The results with Cooper pairing are not very sensitive 
to the core neutrino emissivity. At low accretion rates, the curves are termi- 
nated at the accretion rate where the carbon begins to burn stably (defined as 
the point where depletion time for carbon equals the recurrence time). 



burst sources is believed to lie in the narrow range 0.1-0.3 Ed- 
dington, but there is some uncertainty in these estimates due 
to for example uncertainty in the relation between accretion 
rate and X-ray luminosity, and distance uncertainties. Figure 
^] shows the ignition column depth and the energy per gram 
released in the crust that flows outwards Q/, as a function of 
accretion rate for different models. 

For normal crust neutrons, the constraint on core neu- 
trino emission can be relaxed if the accretion rate is in 
fact larger than inferred from the observed X-ray luminos- 
ity. For example at m » 0.5 «?Edd, standard slow cooling 
from modified URCA reactions explains the observed igni- 
tion columns. However, with Cooper pair cooling in the crust 
(dotted curve in Fig. 1 1 51 . the ignition columns remain well 
above 10 12 g cm" 2 even at accretion rates « niEdd- This im- 
plies that some extra heating of the carbon layer must occur 
that is not included in our model. 

Another important point is that with Xc = 0.2, unstable ig- 
nition requires m > 0.3 niEdd, because at lower rates the car- 
bon burns stably (the curves in FigurelT5lterminate on the left 
where the carbon begins to burn stably). A carbon fraction 
Xc > 0.2 is consistent with the results of our lightcurve fits, 
which gave En > 2 in most cases. 

Predictions for the observable quantities recurrence time 
and superburst energy are shown in Figure ^] We indi- 
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FIG. 16. — For the models shown in Figure fT5l we plot the energy released 
from the surface in the first 6 hours following ignition (as calculated using the 
cooling models of §2), and the recun'ence time. The errorbars indicate the 
estimated recun'ence times and accretion rates for most superburst sources 
and for the near-Eddington accretor GX 17+2. 



cate the observed constraints on recurrence time for super- 
burst sources, and separately for the rapidly accreting source 
GX 17+2. Only the models with low neutrino emissivity in 
the crust and core come close to the observed values at the 
estimated accretion rates. The accretion rate for GX 17+2 is 
quite uncertain; we adopt a value of WEdd for this source. For 
a given ignition depth and carbon fraction, we use our cooling 
models to predict the superburst energy, which we take to be 
the energy released in the surface in the first 6 hours. Note 
that the energies are close to 10 42 ergs in Figure [161 because 
we take Xc = 0.2 in these models, the energy would be sig- 
nificantly smaller if Xc = 0. 1 . The behavior of the superburst 
energy with accretion rate is in general not constraining. For 
long recurrence times > 1 year, the superburst energy satu- 
rates at w 10 42 erg because of the effects of neutrino cooling 
and inwards conduction of heat, as discussed in §2. 

The overall conclusion is that to achieve ignition at column 
depths implied by our fits requires inefficient neutrino cooling 
from the core and the crust, and accretion rates larger than in- 
ferred from the X-ray luminosity. Even with normal neutrons 
and core neutrino emission that is less efficient than modified 
URCA, it is difficult to reproduce the observed superburst re- 
currence times and column depths at accretion rates as low 
as 0.1 /«Edd- Neutrino cooling from the crust due to Cooper 
pair formation results in ignition depths that are too large even 
for accretion rates near Eddington. When this process is in- 
cluded, our models cannot reproduce the observed superburst 
recurrence times. In addition, carbon fractions of > 0.2 are 
required to avoid stable burning of the carbon and achieve un- 
stable ignition at accretion rates of 0.3 mEdd- 



11 



E 



1 1 " 

\ 

\ 

\ 

- \ 


1 1 1 | 1 I I , 1 1 1 | 1 1 1 | I 

— — — e 


~~~~~~~~~ c 
^^C^" " e SF 


\ 

\ 




b 


\ / 


/ / y 

/ / yr 
X / / 






"\/ 






' — \ 

\ /""^ 


a \ 


, 1 , , , 


\ 

1 
1 

. . . i \ . . i . 





8 10 12 14 16 18 

log Column Depth (g cm -2 ) 

FIG. 17. — Temperature profiles for pure helium ignition models at m — 
0.01 njEdd> with a disordered lattice. We show five examples of core neutrino 
emissivity, a to e from Table 2. The dot-dashed curve is model e including 
Cooper pair emission from the crust. The long-dashed curve is the triple 
alpha ignition curve. The dotted line marks a column depth of 10 10 g cm" 2 , 
as inferred for the long burst from 2S 0918-549 observed by in 't Zand et 
al. (2005). 



4. IGNITION MODELS FOR PURE HELIUM BURSTS 

We now consider the constraints that come from pure he- 
lium accretors. Pure helium bursts are interesting because 
they can occur at a wide range of accretion rates. Two sources 
in particular are thought to be accreting pure helium (per- 
haps with a small amount of hydrogen; Podsiadlowski, Rap- 
paport, & Pfahl 2002), with accretion rates different by an 
order of magnitude. The ultracompact binary 4U 1820-30 has 
an orbital period of only 1 1 .4 minutes (Stella, Priedhorsky, 
& White 1987), implying a hydrogen-poor companion. This 
source shows frequent and regular Type I X-ray bursts whose 
properties are consistent with accretion of pure helium at rates 
close to M w 0.1-0.2 Mem as inferred from the X-ray lumi- 
nosity at the time when bursts are seen (Bildsten 1995; Cum- 
ming 2003a). The persistent X-ray source 2S 0918-549 is 
suspected also to be an ultracompact binary because of its 
low optical to X-ray flux ratio, and lack of hydrogen lines in 
its optical spectrum (Juett et al. 2001; Nelemans et al. 2004). 
Recently, a long duration X-ray burst was observed from this 
source whose properties can be explained by accretion of pure 
helium at the observed rate of M « 0.01 MEdd (in 't Zand et 
al. 2005). 

We have calculated ignition conditions for pure helium 
bursts as a function of accretion rate, crust composition, and 
core neutrino emissivity. The calculations follow those for 
superburst ignition in §3, except that the accumulated fuel is 
pure helium, and the nuclear burning rate is given by the triple 
alpha rate from Fushiki & Lamb (1987a). We start integrat- 
ing at a column depth of 10 3 g cm" 3 , and set the temperature 

there proportional to Ft 1 , although the solutions are not very 
sensitive to the outer temperature in most cases. 

We first consider accretion at a rate m = 0.01 mEdd appro- 
priate for 2S 0918-549. Temperature profiles for this accre- 
tion rate are shown in Figure ^1 We show profiles for the 
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FIG. 18. — The outwards heat flux Qi,, ignition column depth y\g n , and 
predicted energy and recurrence times for pure helium flashes, for the same 
models as Figure lTTl In addition, we show a model with Q = 100 in the crust 
for core emission "c" (dot-dashed curve), and a model with Cooper pairs 
included for core emission "e" (dashed line). 



five different core neutrino emissivities in Table 2, and for 
core emissivity "e" but including Cooper pair emission. At 
m w 0.01 /«Edd> Cooper pair emissio n lim its the crust tem- 
perature to T k 2.3 x 10 8 K (see eq. IB6I in Appendix B for 
an analytic estimate). The long burst from 2S 0918-549 (in 
't Zand et al. 2005) had an energy of 10 41 ergs, implying 
an ignition column depth of y w 10 10 g cm" 2 for an energy 
release of w 10 18 erg g" 1 appropriate for helium burning to 
heavy elements. This column depth is also consistent with the 
burst lightcurve, which is well-fitted by a cooling model based 
on this column depth (in 't Zand et al. 2005). FigurefTTlshows 
that ignition at y ~ 10 10 g cm" 2 requires that the core neutrino 
emissivity not be more efficient than modified URCA. En- 
hanced slow cooling (model b), e.g. by Cooper pairing in the 
core, or fast core cooling (model a) lead to ignition at columns 
of > 10 11 gem" 2 . 
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FIG. 19. — Predicted burst energy against recurrence time for pure helium 
flashes. For comparison, we show observed burst properties for 4U 1820-30, 
and 2S 09 1 8-549. The solid curves are for a disordered crust with different 
core neutrino emissivities, the dot-dashed curve is for standard slow cooling 
with Q = 100 in the crust; the dashed curve is for Cooper pairing in the crust 
and core emission "e". 



This conclusion is sensitive to the assumed accretion rate. 
However, models d and e require substantial increases in the 
accretion rate to obtain burst energies « 10 41 ergs, by factors 
of 5-10 above the inferred rate of 0.01 M^&d- Figure[T8lshows 
the variation in the heat flux from the crust Qi,, ignition depth, 
and burst energy and recurrence time with accretion rate. The 
ignition conditions at low accretion rates are not very sensitive 
to the crust composition, and depend mostly on core neutrino 
emission. For comparison with the disordered crust models, a 
crust with Q = 100 and standard modified URCA slow cooling 
in the core is shown as the dashed curve in Figure H~8l 

4U 1820-30 accretes at a rate comparable to the super- 
burst sources, m « 0.2 m^di- The burst properties observed 
in 4U 1820-30 are shown in Figure^] This source under- 
goes periodic variations in accretion rate, with bursts being 
seen in the low state when the accretion rate is a factor of two 
or more below the time-averaged rate. Because the crust tem- 
perature profile is set by the time-averaged rate, we correct for 
this by plotting the burst properties at the time-averaged rate, 
and decreasing the recurrence time by a factor of two. The 
recurrence time is again better explained if the neutrino emis- 
sion is inefficient in the core and crust. Whereas Cooper pair 
neutrino emission does not significantly affect the recurrence 
times at low accretion rates, it does make a difference at the 
higher accretion rates appropriate for 4U 1820-30 because of 
the larger crust temperatures. Our results are consistent with 
the previous analysis of Cumming (2003a), who noted that to 
achieve ignition for pure helium at the observed rates required 
a flux from the crust of Qi, 0.4 MeV per nucleon. An es- 
timate of the residual heat released between bursts was not 
enough to account for this extra flux; inefficient core neutrino 
emission offers a new explanation. Another complication for 
this source is that the accreted material may contain a small 
amount of hydrogen, which can significantly heat the accu- 
mulating layer and shorten the recurrence time (Cumming 
2003a). 



An accurate measurement of recurrence time for the long 
helium burst would impose further constraints. In Figure [191 
we plot burst energy against recurrence time for the different 
models. This plot shows clearly that the long duration bursts 
at low accretion rates are most sensitive to core neutrino emis- 
sion, whereas short bursts at higher accretion rates are most 
sensitive to crust composition. The two systems 4U 1820-30 
and 2S 0918-549 therefore offer complementary constraints 
on the interior model. Unfortunately, the recurrence time for 
the 2S 0918-549 burst is not well constrained (in't Zand et 
al. 2005). In Figure^] we show the observed lower limit of 
1.1 days. 

To summarize, the conclusions for pure helium bursts are 
the same as for superbursts. Enhanced core neutrino emission 
relative to modified URCA leads to larger ignition columns, 
energies, and recurrence times than observed for 2S 0918- 
549. Cooper pairing in the crust leads to much larger recur- 
rence times and energies than observed for 4U 1820-30. The 
constraints again depend on the assumed accretion rate, and 
can be relaxed if the accretion rate is larger than inferred from 
the X-ray luminosity by factors of > 2. 

5. SUMMARY AND DISCUSSION 

In this paper, we have compared models of carbon and he- 
lium ignition on accreting neutron stars to observations of 
long duration X-ray bursts. In particular, we have investi- 
gated the effect of the thermal profile of the crust and core on 
the ignition conditions, and how well the ignition conditions 
reproduce the observed burst properties. We have improved 
on the earlier work of Brown (2004) and Cooper & Narayan 
(2005) by (i) using cooling models of superbursts to predict 
observational properties, and then using these lightcurves to 
provide an independent constraint on ignition depth and ener- 
getics, (ii) including neutrino cooling in the inner crust due to 
Cooper pairing of neutrons, and (iii) considering pure helium 
accretion in addition to carbon. 

5.1. Superbursts 

We applied the cooling models for superbursts calculated 
by Cumming & Macbeth (2004) to observed superburst 
lightcurves. Despite the large uncertainties in distance, we 
find that the energy release and ignition column depths are 
quite well constrained, with the best fitting models giving 
E\-i w 2 and y\2 = 0.5-3. Lower values of £17 give a lumi- 
nosity at early times that is lower than observed, or equiv- 
alently, total superburst energies much smaller than the ob- 
served energies of « 10 42 ergs. An upper limit on En comes 
from the lack of photospheric radius expansion observed in 
most superbursts. Large values of En > 2 lead to extended 
periods of super-Eddington luminosities (durations of min- 
utes and longer), inconsistent with observations. The column 
depth is determined by the rate at which the luminosity falls 
in the tail of the superburst. For example, the range of fit- 
ted values y\2 = 0.5-3 goes from 4U 1636-54 at one end to 
4U 1254-690 at the other. The superburst from 4U 1254-690 
is the longest observed, taking 7 hours to fall to a luminosity 
30% of the peak value, whereas the superburst from4U 1636- 
54 had a much shorter duration, falling to 30% of the peak 
luminosity after roughly 1.5 hours. In our models, large ig- 
nition column depths of « 10 13 g cm" 2 lead to much slower 
decays than observed. 

The fitted column depths are roughly consistent with obser- 
vational constraints on recurrence times. At the moment, only 
4U 1636-54 and GX 17+2 have shown multiple superbursts. 
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FIG. 20. — Ignition column depth for carbon as a function of the flux from 
below heating the accumulating fuel layer. We write the flux as the equivalent 
energy per nucleon at 0.1 of the Eddington accretion rate. The solid curve 
are models with iron as the heavy element, the dashed curve is for a heavy 
composition A = 104, Z = 44. The ignition column is not very sensitive to 
carbon abundance; for this Figure we assume Xq = 0.2. 



For 4U 1636-54 we derived a column depth of yyi « 0.5 by 
fitting the superburst lightcurve. The time between the pre- 
vious superburst and the superburst observed by RXTE/PCA 
was 1.75 years (Kuulkers et al. 2004). Taking this to be the 
recurrence time and combining with the fitted ignition col- 
umn gives an accretion rate m = 0.10 fliEdd, exactly as in- 
ferred from the persistent X-ray luminosity. For GX 17+2, 
the mean recurrence time is 30 days (in 't Zand et al. 2004), 
giving an accreted column depth of 2 x 10 11 g cm" 2 for ac- 
cretion at m = niEdd- The fitted column depth for burst 2 from 
this source is 6 x 10 11 g cm -2 , a factor of 3 larger. The gen- 
eral constraints on superburst recurrence times are 0.4-2 years 
(in't Zand et al. 2003), which gives column depths y» 1- 
5 x 10 11 gem" 2 (m/0.1 m Edd ). 

How do the fitted column depths compare to carbon igni- 
tion models for superbursts? The most striking result is that 
to achieve ignition at the inferred column depths for accre- 
tion rates thought to be appropriate for these sources, 0.1- 
0.3 Mndd, requires adjusting each parameter to maximize the 
flux emerging from the crust. To illustrate this, Figure |20] 
shows the ignition column depth as a function of the base flux. 
We find that the flux required for ignition 8 at y = 10 12 g cm" 2 
is Qt ~ 0.25 (ffi/0.3 /«Edd) _1 - FigurefTBIshows that this value 
of flux requires that (i) crust cooling by Cooper pairs is not 
active, (ii) the core neutrino emission is significantly reduced 
below modified URCA, and (iii) the crust conductivity should 
be "disordered" so that it has a low thermal conductivity. We 
now discuss the likelihood of satsifying each of these require- 
ments. 

Most important is the effect of Cooper pairs in the crust. 

8 We give the result here for a heavy element of 56 Fe. The argument for a 
heavier composition is similar because as Brown (2004) pointed out, the flux 
emerging from the crust is smaller, compensating for the fact that a smaller 
flux is needed for ignition at a particular column depth. Nonetheless, a heav- 
ier composition does give a smaller ignition column depth, but the crust and 
core properties are more important parameters. 
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FIG. 21. — Maximum outwards flux and minimum column depth for un- 
stable carbon ignition for different carbon fractions. At a given accretion 
rate, a larger flux than indicated in the upper panel results in stable rather 
than unstable burning of carbon. This translates into the minimum possible 
column depth for unstable carbon ignition shown in the lower panel. The 
dotted lines indicate the range of accretion rates for most superburst sources, 
0.1-0.3 m EM . 



The emissivity due to this process is large wherever the tem- 
perature of the crust is close to the critical temperature for 
superfluidity T c , because this allows the efficient formation 
and breaking of Cooper pairs (Yakovlev et al. 1999). This 
happens near the base of the crust, but most importantly close 
to neutron drip where most of the nuclear energy release in 
the crust occurs. The effect is to limit the temperature of 
the crust to be < 4.4 x 10 8 K (m/m Edd y/ 7 (eq. ESJ). In Ap- 
pendix C, we discuss the uncertainties in this cooling mecha- 
nism. The emissivity depends on the T c profile with density, 
however all models that we have tried give substantial cool- 
ing rates from this process, resulting in ignition columns for 
superbursts >4x 10 12 g cm" 2 for accretion rates less than the 
Eddington rate. Some extra heating of the accumulating car- 
bon layer is needed to explain observed superburst properties 
if Cooper pair cooling is active. 

Suppression of the core neutrino emissivity below the mod- 
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FIG. 22. — Ignition column depth for pure helium as a function of the flux 
from below heating the accumulating fuel layer. We write the flux as the 
equivalent energy per nucleon at 0.1 of the Eddington accretion rate. 



ified URCA rate will occur if either the protons or neutrons 
are superfluid in the core. However layers whose tempera- 
ture T < T c will copiously emit neutrinos through the Cooper 
pair process and hence reduction of the core neutrino luminos- 
ity requires critical temperatures for neutrons and/or protons 
wich satisfy either T c < T or T c ^ T everywhere in the core. 
If both neutrons and protons are paired in the whole core with 
r c 's well above 10 9 K everywhere, then only minor neutrino 
emission processes may be active and the model (e) of Table 
12 would apply. Nevertheless, in this case the star's specific 
heat would be very reduced leading to rapid cooling of iso- 
lated cooling neutron stars of age around 10 5 years, in conflict 
with observations (see, e.g., Gusakov et al. 2004). Yakovlev 
& Pethick (2004) argue that the best calculations of transi- 
tion temperatures imply that the protons are superfluid, but 
neutrons not. In the case of proton superconductivity, and as- 
suming T c ^> 10 9 K, the core neutrino emissivity would be 
set by neutron-neutron Bremsstrahlung, roughly an order of 
magnitude below the modified URCA rate, corresponding to 
model (d) of Table |2 However, all calculations of proton 
gaps show that T c decreases and eventually vanishes at rel- 
atively low density (see, e.g., Figure 9 in Page et al. 2004) 
and hence strong proton superconductivity, to avoid the pro- 
ton 'So Cooper pair neutrinos, probably requires a very low 
mass star. Moreover, considering that all superbursters are 
heavily accreting, and probably also very old, some sources 
may be massive stars in which case, due to the resulting high 
central density, numerous fast neutrino processes of the direct 
URCA family are possible (see, e.g., Page et al. 2005 for a 
review). 

Finally, the crust conductivity is expected to be low in ac- 
creting neutron stars, because hydrogen and helium burning 
produces a complex mixture of heavy elements (Schatz et 
al. 1999). Our results suggest a completely disordered crust 
with a thermal conductivity essentially equivalent to that of 
the liquid state. Even a very impure crust with impurity pa- 
rameter Q = 100 does not fit the data as well as a completely 
disordered crust. Interestingly, recent work by Jones (2001, 
2004a, 2004b) concludes that the same is likely to be true for 
the original crust before accretion starts. 



It is very important to stress that conclusions about the in- 
terior thermal properties are sensitive to the choice of local 
accretion rate. The constraints on core emissivity can be re- 
laxed by an increase in accretion rate by factors of 2 or 3. This 
is a possibility since the relation between X-ray luminosity 
and accretion rate is uncertain, and the accreted material may 
cover only part of the stellar surface (Bildsten 2000). How- 
ever, we stress that if Cooper pair cooling operates, we cannot 
reproduce observed superburst properties even for accretion 
rates near the Eddington rate. 

The carbon ignition models also show that a large carbon 
fraction Xc > 0.2 is needed if conditions for the thermal in- 
stability are to be achieved before the carbon stably burns 
away. This is illustrated in Figure |^ which shows the max- 
imum value of Qb that allows stable burning as a function 
of m and X c - For larger base fluxes, the carbon burns sta- 
bly. This limit on base flux translates into a lower limit to 
the ignition column depth at particular accretion rate, shown 
in the lower panel of Figure |^ Ignition at column depths 
near 10 12 g cm" 2 at the observed accretion rates for superburst 
sources requires Xc > 0.2. This conclusion is consistent with 
the lightcurve fits, which imply En ks 2, if carbon burning 
releases m 10 I8 Xc erg g" 1 as expected for carbon burning to 
iron group. There is no need for additional energy release, for 
example from photodisintegration of heavy elements (Schatz, 
Bildsten, & Cumming 2003a). The large inferred carbon frac- 
tion is an important constraint on models of rp-process hydro- 
gen and helium burning. Current models suggest that stable 
burning of hydrogen and helium most likely plays a role in 
producing the fuel (in 't Zand et al. 2003; Schatz et al. 2003b), 
although further work is needed. 

5.2. Pure helium bursts 

Pure helium bursts are interesting because they probe lower 
ignition masses and densities than superbursts, and can occur 
at a wider range of accretion rates. The likely ultracompact bi- 
nary 2S 0918-549 is a persistent source as far as is known, and 
accretes at a rate « 0.01 Mem- By considering the energetics 
and by fitting the burst lightcurve with the cooling models of 
CM04 extended to low column depths, in 't Zand et al. (2005) 
showed that this burst is consistent with pure helium ignition 
at y « 10 10 g cm" 2 . They also showed that pure helium accre- 
tion at ss 0.01 MEdd gives ignition at 10 10 g cm" 2 if most of the 
heat released in the crust flows outwards. Figure|22]shows the 
ignition column depth for pure helium accretion as a function 
of base flux. At O.O^Edd, Qb ~ 1 MeV per nucleon is re- 
quired for ignition at y f» 10 10 g cm" 2 . Our ignition models in 
this paper show that this requires that the core neutrino emis- 
sivity not be enhanced over modified URCA. Either a slow 
cooling rate enhanced by a factor of 30, or fast cooling in the 
core give ignition column depths > 10 11 g cm" 2 at this accre- 
tion rate. An increase in accretion rate of factors of > 5 over 
the assumed value is required to bring the enhanced cooling 
models into agreement. At low accretion rates, the ignition 
depth is most sensitive to core temperature, and is not very 
sensitive to crust properties. In particular, the crust neutrino 
emission plays a smaller role than for superbursts because of 
the lower crust temperatures. 

Pure helium bursts are also observed at a similar accre- 
tion rate to superburst sources, from the ultracompact binary 
4U 1820-30. In this case, we showed that the observed burst 
properties are again best fit by models which maximize the 
outwards flux from the crust. Models with Cooper pair emis- 
sion from the crust give ignition depths, recurrence times and 
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energies that are too large by a factor of 5. This is consistent 
with the previous conclusions of Bildsten (1995) and Cum- 
ming (2003), who did not consider the crust or core physics, 
but noted that the helium layer must be quite hot to achieve ig- 
nition at the depth inferred from observations. If the accretion 
rate is larger by a factor of > 2 than inferred from the X-ray 
luminosity, these constraints are relaxed. The uncertainty as- 
sociated with the accretion rate can be bypassed by studying 
the burst energy as a function of recurrence time, as shown in 
Figure^] however, current constraints are limited. Unlike at 
low accretion rates, the ignition conditions mainly depend on 
crust properties at high rates, giving a complementary view of 
the interior. 

5.3. Conclusions and Future Work 

The observational progress on long Type I X-ray bursts 
has opened up a new probe of accreting neutron star in- 
teriors, complementary to studies of isolated neutron stars 
(e.g. Yakovlev & Pethick 2004, Page et al. 2005) and accret- 
ing neutron stars in quiescence (Brown, Bildsten, & Rutledge 
1998; Colpi et al. 2001; Rutledge et al. 2002; Wijnands et 
al. 2002; Yakovlev et al. 2004). We find in this paper that the 
long Type I burst from 2S 0918-549, a pure helium accretor at 
0.01 MEdd, is best fit by models with core neutrino emissivity 
equivalent to modified URCA or smaller. For superbursts, and 
pure helium bursts from 4U 1820-30, which occur at higher 
accretion rates > 0. 1 the ignition models limit both core 
and crust neutrino emission. In particular, neutrino cooling by 
Cooper pairing of neutrons in the crust leads to superburst ig- 
nition column depths that are too large. Either the Cooper 
pairing emissivity is much less than current calculations sug- 
gest 9 (see Appendix C for an assessment of the uncertainties 
in this process), or an additional heating source not included 
in current superburst ignition models is required in the accu- 
mulating fuel layer. 

Our models can be improved in several respects. First, 
our ignition models use the one-zone criterion of Fujimoto, 
Hanawa, & Miyaji (1981) and Bildsten (1998) to estimate 
the ignition column depth. Although this technique compares 
well with numerical simulations and normal mode analyses 
(e.g. Woosley et al. 2003; Cooper & Narayan 2005), time- 
dependent calculations of ignition should be carried out to cal- 
culate ignition conditions. These calculations are in progress 
(Halpin & Cumming 2005, in preparation). Our cooling mod- 

9 An alternative explanation is that these sources are not neutron stars but 
rather "strange stars" (e.g. Alcock, Farhi, & Olinto 1986). Strange stars do 
not have an inner crust which would naturally explain the lack of emission 



els for burst lightcurves assume instantaneous burning of the 
fuel, and cannot address the physics of the rise. Further stud- 
ies of the superburst rise are needed. An additional motiva- 
tion for this is to understand the observed precursors to super- 
bursts, possibly due to triggering of an overlying H/He layer 
by the carbon runaway. 

The thermal models of the interior used here assume steady 
state accretion. In fact, superbursts have been observed from 
transient systems, in which case the core temperature will be 
lower than assumed in our models. Observations of quiescent 
cooling in KS 1731-260 (Rutledge et al. 2002; Wijnands et 
al. 2002), interpreted as cooling of the crust, imply a cold core 
and high crust conductivity, exactly opposite to the conclu- 
sions from superburst ignition calculations, as emphasized by 
Brown (2004). The time-dependent calculations of Rutledge 
et al. (2002) for this source indicate that the crust tempera- 
ture reaches a maximum value of « 2.5 x 10 8 K, lower than 
the temperature in our steady-state models. Cooper pair cool- 
ing in the crust was not included in the models of Rutledge 
et al. (2002), but is probably not important at these low tem- 
peratures. The lower crust temperatures may make achieving 
ignition at the inferred column depth for KS 1731-260 diffi- 
cult, implying that an extra heating mechanism is required in 
superburst models. On the other hand, this may be related to 
the larger ignition column depth for KS 1731-260 compared 
to sources such as 4U 1636-54 for example. Further work 
on the time-dependent thermal profile in transiently accret- 
ing sources and the consequences for superburst ignition is 
needed. 
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due to Cooper pairing of neutrons. This scenario is explored in Page & Cum- 
ming 2005). 



APPENDIX 

A. THE EARLY PHASE OF THE SUPERBURST LIGHTCURVE 

When fitting our time-dependent cooling models to observed superburst lightcurves, the inferred value of Eyj depends on the 
early part of the cooling curve, for times shorter than the thermal time of the fuel layer. Therefore it is important to understand 
the physics of this phase of the lightcurve. In this Appendix, we describe a simple steady-state model of the early cooling which 
highlights the physics, and gives us confidence in our numerical results. 

CM04 showed that after the fuel burns, the early phase of the superburst lightcurve is set by a cooling wave which propagates 
inwards from the surface. At column depth y, there is a characteristic thermal timescale fthermOO ~ H 2 /D, where H is the pressure 
scale height, and D the thermal diffusivity. The thermal timescale grows with depth. After time f , the cooling wave has penetrated 
to a depth at which t « fthermCyi)- At lower column depths y < y/,, the atmosphere has a constant heat flux with depth; at higher 
column depths y > yt,, the temperature profile has not evolved significantly from its initial state. 

This picture suggests a simple model of the early phase of the lightcurve. We first make an analytic estimate for constant opacity 
in the layer, and then present numerical calculations integrating the true opacity profile. For constant opacity, the temperature 
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profile in the constant flux region is given by integrating F = (AacT /3n)(dT /dy) from the surface. The radiative zero solution is 

1/4 



r-WxufKJSra^jjgg^ , (A!) 

where we choose a typical value for the opacity. Equation (lAli gives the temperature profile in the constant flux region. However, 
at the base of the constant flux layer, the temperature remains close to its initial value, which is set by the nuclear energy release 

at the beginning of the flash. CB01 and CM04 calculated this to be 7) = 4.0 x 10 9 K yjf £ , 1 1 7 /2 (g 14 /2.45) 1/s . Combining these 
two conditions by using the radiative zero solution to set T(y/,) = 7] gives the flux through the surface when the cooling wave has 
penetrated to depth y^,, 



0.02 cm 2 g" 1 / V2.45 

All that remains is to relate the depth of the constant flux region to the time after ignition. The thermal time is f t herm ~ H 2 /D, 
where D = K/pcp is the thermal diffusivity, giving 

_ pcpH 2 _ -icpKy 2 

ftherm - R - ^3 , (A3) 



or taking k to be constant, 

f therm = 11 hours ( ?T7 ^W) (^) _1/4 - (A4) 
En \0.02cm z g 1 J V2.45/ 

Now setting t = ftherm and rewriting equation \A2l in terms of ftherm using equation \AA\ specifies the flux evolution with time. We 
find 

F^2.2 (l)" /3 £ yf \ r y 2/3 (^f 12 , (A5) 
Vhrs/ 17 V0.02cm 2 g-V V2.45/ ' 

which compares favorably to the empirical fit to the time-dependent simulations from CM04, F24 ~ 2 (f/hr) -0 - 2 ^ 7 ^ 4 . 

In fact, the opacity varies through the layer, typically by a factor of 5. We therefore repeat the same argument, but integrating 
numerically. For a given column depth yj, we find the flux F required so that a constant flux atmosphere has a base temperature 
which matches the initial temperature at that depth (set by the initial energy release). We then calculate the thermal time at co lumn 
depth yt,, which tells us the time at which the cooling wave reaches yj or when this value of surface flux applies. Figure lA23l 
shows the resulting lightcurves, compared with our time-dependent cooling models. The agreement is good, giving us confidence 
in the results of our numerical models. 

It is important to note that the radiative op acity c an be important even close to the base of the burning layer 10 , because of 
the high temperatures reached initially. Figure lA24l shows the opacity profile in the layer immediately after the fuel burns. The 
radiative opacity controls heat transport at lower densities, and has contributions from electron scattering and free-free absorption. 
At high density, the growing free-free opacity shuts off the radiative heat transport, and electron conduction takes over. As the 
layer cools, the transition from radiative to conductive cooling happens at lower and lower densities, since K con d oc T 2 , and 
Kff oc r~ 7 ' 2 , moving the intersection point to lower density. A similar transition was noted in models of steady H/He burning by 
Schatz et al. (1999), but at lower column depths because of the lower temperatures (5-10 x 10 s K, in that case). For the free-free 
Gaunt factor, we use the fit of Schatz et al. (1999) to the calculation of Itoh et al. (1991). We have checked that the physical 
conditions in our models are within the regime of applicability of these calculations. 

B. ANALYTIC MODELS OF THE CRUST TEMPERATURE PROFILE 

Figure [D^] shows that the flux entering the accreted layer (and therefore the burst ignition conditions) mainly depends on the 
core neutrino emissivity at low accretion rates, and the crust properties at high accretion rates. Brown (2000) computed the 
thermal structure of rapidly accreting neutron stars, and pointed out that the thermal structure of the crust could be understood 
using a simple analytic model. We take a similar approach here to understand the variation in the temperature profile as a function 
of accretion rate. 



Without neutrino cooling in the crust 

We first ignore neutrino emission from the crust, and assume that the heat released in the crust Q nuc w 1.4 MeV per baryon is 
either radiated from the surface or conducted into the core and radiated as neutrinos. For accretion rates m > 0.1 m^dd, most of 

10 The plasma frequency u) p in the burning layer is quite large, and one might therefore expect a suppression of radiative heat transport caused by the inability of 
photons with u) < (dp to propagate (we thank D. Yakovlev for raising this issue). This is expected for temperatures T <tiT p , where T„ = Huip/ks ~ 10 10 K (pgY e ) 1 / 2 . 
We have calculated the suppression factors for our models based on the modified Rossland mean opacity integrals given by Aharony & Opher (1979) and van 
Horn (1992), and find that this effect is not important, changing the calculated cooling curves by of order 1%. The smallness of the effect is firstly due to the 
fact that T is always close to or greater than T p when radiative transport dominates electron conduction (at lower temperatures, T p /T is larger, but radiative heat 
transport is no longer important relative to conduction), and secondly because free-free absorption dominates the opacity in the region most likely to be affected 
(its l/u> 3 frequency dependence leads to a preference for higher energy photons than electron scattering in the Rosseland mean integral, and so less suppression 
at a given temperature; Aharony & Opher 1979). 
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FIG. A23. — Flux as a function of time for £17 = 1 (lower curves) and £17 = 2 (upper curves). The solid curves show the simple steady-state model described 
here; the dotted curves are the results of time-dependent simulations (with y l2 = 1). 
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FIG. A24. — Opacity (continuous line) as a function of depth immediately after burning for £17 = 2. We also show the separate contributions to radiative 
opacity from electron scattering (es; dotted curve) and free-free absorption (if; short-dashed curve), and the contribution from electron conduction (cond; long- 
dashed curve). At low densities, electron scattering dominates the opacity; at high densities, free-free absorption blocks radiative transport of heat, and electron 
conduction takes over. 



the energy enters the core. In this case, the peak temperature in the crust is whatever it needs to be to conduct the heat into the 
core, given the thermal conductivity. If electron-ion collisions set the conductivity, as for a disordered crust, then 



ii 2 k\Tn e \ I 3nh 



s3 



K \ 3m, ) \4e*m*Z j (B1) 

(e.g. Yakovlev & Urpin 1980), where we assume a single species of ions with charge Z, m + = E F /c 2 is the electron effective 
mass, and the second term is the inverse of the electron-ion collision frequency. We have set the Coulomb logarithm to unity, a 
reasonable approximation for our purposes. 

We now solve F = pKdT / dy in the inner crust to find the temperature profile. In the inner crust, the pressure is mostly from 
degenerate non-relativistic neutrons, giving P oc (Y n p) 5 / 3 . Using this to integrate, we find 

where we write the fraction of the energy released in the crust that flows inwards as To write this expression, we have assumed 
that the core temperature is much smaller than the maximum crust temperature, and that the density at the crust/core boundary 
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is much greater than the density at the location of the temperature maximum (these approximations are good enough for our 
purposes). We have also taken Y„ = 0.8, Y e = 0.05, and Z = 22 in the inner crust. The accretion rate enters this formula because 
the inwards flux is F = f m Qmx<jh. 

Similarly, we can integrate into the crust from low densities. In the outer crust, degenerate relativistic electrons set the pressure. 
In this case, the temperature profile is given by 



r 8 « 7.9 ^ 



0.1/ I 12 J lm Ed J 



HP/Po) 



1/2 



(B3) 
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(see also CB01), where Pq is the outer pressure at which we start the integration, and we have assumed that the temperature there 
is small. 

Now matching the outer and inner temperatures, and assuming that / ut + /m = 1 and / ou t < 1, we find 

Q fe «e nu c/out~0.3MeVpernucleon (j^j ■ ( B4 ) 

The prefactor in this estimate is approximate. Brown (2000) gives a mo re careful analysis using more accurate power laws for 
the equation of state in the outer and inner crust. Howeve r, eq uation JB4> nicely reproduces three properties of Qb that we find in 
the numerical c alcul ations at m > 0.1 mEdd- First, Fi gures \l5\ and[T8l show that at these accretion rates, Qb is not very sensitive to 
m. In equation JB4L m drops out altogether because the temperature profiles in the inner and outer crust both scale in the same 
way with m. Second, Qb is insensitive to the core neutrino emission as long as the maximum temperature in the crust is much 
greater than the core temperature. Thirdly, equation flJ34b shows that Qb is inversely proportional to Z 2 /A. This was pointed out 
by Brown (2004) (the actual scaling is Qb oc (Z 2 /A)" 0,6 ), and is the reason why the ignition conditions are much less sensitive to 
composition than found by CB01 (who varied Z 2 /A but kept Qb constant) 11 . 

At low accretion rates m < 0. 1 m^n, a large fraction of the heat released in the crust flows out rather than in, / out ~ 1 , as shown 
in Figure^] The inner crust is almost isothermal in this case (see Fig. 1171 ). and the outwards flux is determined by the core 
temperature, which in turn is set by the core neutrino emission. This makes Qb more sensitive to the core neutrino emission, and 
less sensitive to crust properties, than at higher accretion rates. 

With neutrino cooling in the crust 

We have ignored neutrino emission from the crust. However, as we have seen in the numerical calculations neutrino emission 
can have a dramatic effect on the temperature profile. One way to estimate when crust neutrino emission matters is to ask at what 
temperature the crust neutrino emission balances the energy release in the crust. For crust Bremsstrahlung, the neutrino emissivity 
is Q v w 0.3 erg g" 1 s" 1 TgL(l -Y„)(Z 2 /A), where L is the Coulomb logarithm, and is insensitive to depth (see Fig.^Jand Haensel 
et al. 1996); we estimate a typical value Q v w 3 x 10 12 erg cm" 3 s" 1 T%, and a total crust luminosity L„a3 x 10 30 erg s" 1 T%. 
The total nuclear energy release in the crust is L crast « 10 36 erg s" 1 (m/'"Edd)- Therefore L v = L aust when 

T 8 «8 (m/m E dd) 1/6 , (Bremsstrahlung) (B5) 

or T% « 4 for 0.01 mEdd, and 7s w 6.5 for 0.3 »?Edd- This is in good agreement with Figures and [n] where neutrino emission 
plays an important role only for the hottest models in these figures. For Cooper pair neutrino emission in the crust, Q v w 
7 x 10 2l) erg cm" 3 s" 1 r g 7 (Icf /fm" 1 ) (Yakovlev et al. 1999). Using a height of 0.5 km for the emitting region gives L u « 3 x 
10 31 erg s" 1 r 8 7 . Again setting L v =L crust , we find that neutrino emission will dominate the energy loss at the lower temperature 

T 8 w 4.4 (m/m Ed d) 1/7 - (Cooper pairs) (B6) 

This is in good agreement with Figure ^] Ignition of carbon at column depths < 10 12 g cm" 2 requires temperatures close to 
6 x 10 8 K at the ignition point (see ignition curve in Fig. II 1>. This is close to the maximum temperature for m = 0.3 ifiEdd without 
Cooper pairs, explaining why we need to go to higher accretion rates to get ignition at y n < 1 . With Cooper pairs, the temperature 
at the ignition point is constrained to be w 4 x 10 K, giving ignition columns of y^ ~ 10. 

Limits on the outwards flux and ignition column depth 

The fact that the crust temperature is limited by neutrino emission from the crust g ives an uppe r limit to the o utwa rds flux 
at a given accretion rate. Substituting the limiting temperature from either equation dB5b or dB6i into equation JB3I for the 
temperature profile of the outer crust, gives 
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(Bremsstrahlung) (B7) 
-i 

(Cooper pairs) (B8) 



where Pi/Pq is the ratio of inner to outer pressures. These limits on the flux translate into limits on the ignition column depth, by 
comparison with Figure|20| The limit on Qb from Cooper pairs severely constrains the carbon ignition depth. If Qbim/riiEiA) < 
0.04, then y; gn > 5 x 10 12 g cm" 2 for carbon. Cooper pair emission cools the crust so efficiently that its temperature is too low to 
achieve the required outwards heat flux for ignition at y w 10 12 g cm" 2 inferred from observations of superbursts. 

Note that an implicit assumption in this argument is that the composition of the inner crust does not change as Z 2 /A in the outer crust changes. In fact, this 
assumption is borne out by the calculations of Haensel & Zdunik (2003). 
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C. NEUTRON 'So PAIRING AND NEUTRINO EMISSION BY THE FORMATION OF COOPER PAIRS 

Neutron 1 Sq pairing 

Given the importance of the neutrino energy losses due to the Cooper pair formation process found in this work, we describe 
in this Appendix in some detail the nature of the underlying physics and robustness of the ingredients needed for our models. 
First of all, the existence of a neutron drip regime is beyond doubt and results from the finite depth (~ 50 MeV) and finite width 
(~ nucleus diameter) of the nuclear potential, i.e., this potential can only accomodate a finite, and small, number of bound states. 
With increasing density, the enormous Fermi energy of the electrons leads to neutronization and the resulting large number of 
neutrons cannot be accomodated within the bound levels, i.e., neutrons have to drip at high enough density 12 . These dripped 
neutrons form a degenerate Fermi liquid and, according to the Cooper theorem (Cooper 1956), they will unavoidably pair and 
become a superfluid if there is any attractive interaction between them (immediately afer the development of the BCS theory 
it was proposed that nucleons in nuclei must pair due to this mechanism; Bohr, Mottelson, & Pine 1958). The important, and 
delicate, issue is the value of the pairing critical temperature T c which depends very sensitively on the strength of the pairing 
interaction. At the density range relevant for a neutron star crust the dominant attractive interaction between the dripped neutrons 
is in the spin-angular momentum channel 'So, which, fortunately is very well understood in vacuum. 

In a medium, many-body effects alter the interaction and much effort has been dedicated to studying them (see, e.g., Lombardo 
& Schulze 2001 for a review). At the densities corresponding to neutron star cores, there are still very large uncertainties in the 
size and density extent of the gaps, but in the low density regime of the crust the situation is much better. Many-body techniques 
used to calculate the size of the neutron 'So gap have increased in sophistication with time and results obtained in the last fifteen 
years show a clear convergence. The most important effect turns out to be the polarization of the med ium which, in some sense, 
screens the interaction between neutrons and results in a reduction of the gap. We plot in Figure lC25l results of the most reliable 
calculations (as well as one example of a calculation which explicitly did not i nclu de the effects of medium polarization for 
illustration of the importance of this effect). Numerical values are listed in Table lC3l as well as reference to the relevant works. 
Both the table and the figure illustrate the present uncertainties on T c . A good indication of the convergence of the models is 
that the most recent calculations of Schwenk, Friman, & Brown (2003) incorporated significant new improvements and obtained 
values for T c very close to previous calculations. Given the range of temperatures in accreting neutron star crusts, 1 - 8 x 10 8 K, 
and considering that the Cooper pair neutrino process becomes negligible when T < 0.2 T c (see below), the important region is at 
neutron Fermi momenta smaller than 0.5 fur 1 , i.e., densities below 10 13 g/cm 3 (Fig lC25t where all T c curves we show are very 
close to each other. At higher densities the differences become significant but, fortunately, they do not seriously affect our results. 
For very low neutron densities one can use the following analytical solution (Schwenk 2004), which incorporates the medium 
polarization effect to calculate the zero temperature energy gap (see next paragraph), 

A(0) = (2/e) 7 / 3 e F exp(7r/2k F a n „) (CI) 

where a nn = -18.5 fm is the neutron-neutron scattering length and e F = k F /2m„. The above formula is formally valid only when 
|&F<2nn| <C 1 but gives reasonable results even at k F ~ 0.1 fm"'. At all densities it is customary to calculate the critical temperature 
T c using the BCS result 

A(0) = 1.76T C . (C2) 
Even though this is not necessarily correct for calculations beyond the BCS approximation, it gives a credible value for T c . 

Neutrino emission by Cooper pair formation 

The essence of pairing is the formation of a condensate of Cooper pairs, with the result that the energy of single-particle 
excitation changes from 



n 2 k 2 nhj we 

m= 2M n -2M n =VFh(k - kF) Wlth VF = hdk 



(C3) 

k=k F 



in absence of pairing to 



ep(Jt) = V(v F n(k-k F )) 2 +A(T) 2 . (C4) 

in presence of pairing. The breaking of a Cooper pair therefore requires an energy of at least 2A(T) while its formation generates 
that much energy. At zero temperature there are no single-particle excitations, but at finite temperature T < T c the neutron fluid 
consists of two components, paired neutrons and excited ones, which are in thermal equilibrium and there is constant formation 
and breaking of Cooper pairs. The energy liberated by the formation of a Cooper pair can be emitted in a neutrino-antineutrino 
pair (Flowers, Ruderman, & Sutherland 1976; Voskresensky & Senatorov 1987). For neutron 'So pairing, the emissivity is 
(Yakovlev, Kaminker, & Levenfish 1999) 



Ql s ° = 1.5 x 10 
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hkf 



M„c 



TjF(T/T c ) ergcnfV . (C5) 



The control function 13 F is plotted in Figure IC261 The shape of F reflects the fact that when T ~ T c only a small fraction of 
neutrons are in the paired state and the gap is small, resulting in a diminished emission of neutrinos with little energy, while when 

12 The only exception to this can occur in a "strange star" made of self-bound strange quark matter, in which normal baryonic matter can exist only below the 
neutron drip density in the form of a thin outer crust (Alcock, Farhi, & Olinto 1986). 

13 Our normalized function F is related to the F of Yakovlev et al. (1999) by F = F/F nm with F max = 4.313. 
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TABLE C3 

Critical temperature T c , vs neutron Fermi momentum k F , for five reliable calculations 



k F [fin" 1 
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AWP II & AWP III: Ainsworth, Wambach, & Pines (1989), WAP: Wambach, Ainsworth, & Pines (1993), CCDK: Chen et al. (1993), SFB: Schwenk, Friman, & 
Brown (2003). 

t For values of kf this small, T c cannot be found from the above references, except in the "SFB" case: using Eq. lCll one obtains T c = 1.9 X 10 8 K (Schwenk 
2004). 
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FIG. C25. — Modern results for the neutron 'So pairing critical temperature T c as a function of the neutron Fermi momentum kf, from Table IC3l The curve 
labelled "CCKS" from an older calculation (Chen et al. 1986), in contradistinction to the other ones does not include medium polarization effects and is shown 
for illustration. The upper scale shows the corresponding densities for a crust chemical composition as used in the present work 



T <C T c it becomes impossible to break pairs and hence pair formation does not occur anymore, resulting in a strong suppression 
of the emissivity approximately by a Boltzmann-like factor exp-A(0)/A.flr. 

This neutrino emission process is extremely efficient but not "exotic" in any way and has to be considered in any cooling 
neutron star model, even minimal (Page et al. 2004). Finally, we mention that a similar process occurs in the laboratory when 
free electrons are injected into a superconductor, although the energy released is emitted predominantly in phonons (Schrieffer & 
Ginsberg 1962) instead of neutrinos (!), leading to a pair formation rate in agreement with experiments (see, e.g., Carr et al. 2000 
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for recent results). 
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